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' ' ABSTRACT 

fvj GRS 1915+105 is a very peculiar black hole binary that exhibits accretion-related 

►^ states that are not observed in any other stellar-mass black hole system. One of these 

,-_l states, however - referred to as the plateau state - may be related to the canonical hard 

Q^ state of black hole X-ray binaries. Both the plateau and hard state are associated with 

T-H steady, relatively lower X-ray emission and flat/inverted radio emission, that is some- 

CN times resolved into compact, self-absorbed jets. However, while generally black hole 

f-^ binaries quench their jets when the luminosity becomes too high, GRS 1915-1-105 seems 

f~^ to sustain them despite the fact that it accretes at near- or super-Eddington rates. In 

f^ order to investigate the relationship between the plateau and the hard state, we fit 

T-H two multi-wavelength observations using a steady-state outflow-dominated model, de- 

IL^ veloped for hard state black hole binaries. The data sets consist of quasi-simultaneous 

observations in radio, near-infrared and X-ray bands. Interestingly, we find both sig- 
nificant differences between the two plateau states, as well as between the best-fit 
^ model parameters and those representative of the hard state. We discuss our interpre- 

tation of these results, and the possible implications for GRS 1915-|-105's relationship 
to canonical black hole candidates. 

Key words: Black hole physics - accretion, accretion discs - radiation mechanisms: 
general - X-rays: binaries - galaxies: active - galaxies: jets 



1 INTRODUCTION hole binary (BHB) discovered to ex hibit superliuninal ve- 

^, . ^T,ri,„ir-,„^ ,. , ,^A locities in its radio emitting-ejecta (IMirabel & Rodriguezl 

ine microquasar GKb 1915-1-105 was discovered on 15 Au- ■ ■ - - ' : ^—- ' 

gust 1992, by the WATCH all-sky monitor, aboard the Rus- 



X 



cd 



19941. The obvious parallels to the jets in Active Galactic 



sian GRANAT satellite ( Castro- Tirado, Brandt, & Lund 



Nuclei (AGN) led to this source being classified as a "mi- 



sian i^riAiNAi sateuite v^astro-iiraao, oranax;, n ijuna „ ir-p — 1_ i c — ?> — T~- i mno k a -l. ^ u 

I r=i -— — ; ,' II .1 — - — : — ^ T-^ 1 croquasar JMirabel & Rodriguezll998p. subsequent obser- 

1992 Castro-Tirado et al. 1994h. It is a hard X-ray tran- /^ -.^nW— ;^ 7 r^ — J— r— '„ - v^ ^. - 

1 1 ■ ■ IL r ... , „ vations with instruments onboard the Hossi X-ray I immg 

sient located m the constellation or Aquila, at I = 45.37 , 



Explorer (RXTE) have revealed a richness in variability, dis- 
= —0.22 , and was the hrst stellar mass accreting black ,. . , . /^no irvir , ir,r r ii i -mr-n 

" tmguishmg GRb 19154-105 irom every other known BHB, 

over which astronomers are still puzzling to this day. 
* In Memoriam The longer term X-ray variability exhibits "dipping" 



P.C.N, van Oers et al. 



that is thought to be associated with the recession and re- 
generation of the inner parts of the accretion disc, perhaps 
caused by the onset of thermal-viscous instabilities (iBelloni 



et al.|[T997b|al. Other models have interpreted the spectral 



changes as resulting from the disappearance of the corona 
| Chaty|[l998l [Rodriguez et al.||2008a|b| ), or from the dis- 



sipation of magnetic energy via magneto-hydrodynamical 
plasma processes (e.g. 'Tag ger et al.|2004 |. 

Beyond the dipping behaviour, Belloni et al. ( 2000 1 were 



able to classify all variability patterns stretching over more 
than a year into only twelve classes ( Klein- Wolt et al.|2002 



Hannikainen et al. 2005 later identified two more classes), 
based on colour-colour diagrams and light curves. Ten of 
these twelve classes can be understood as the interplay of 
two or three of three basic states, designated as states A, 
B and C. The remaining two classes, (j) and x, do not show 
state transitions and appear exclusively within states A and 
C respectively. State A displays the highest flux and the 
softest spectrum while state C displays the lowest flux and 
the hardest spectrum. Although state B never lasts for more 
than a few hundred seconds, </> and x episodes can persist 
for days or short intervals of < 100 s. 

Aside from the existence of so many distinct accretion 
states, GRS 1915-1-105 appears similar to other BHBs. Thus 
there has been much discussion (e.g. |Reig, Belloni, fc van der| 
|Klis '2003) about the extent to which any of the states found 
in GRS 19154-105 correspond to the "canonical" states (see 
e.g. [Homan fc Belloni|[2505l [Remillard fc McClintock|[2bo6l 
for deflnitions) found in the average BHB. BHBs generally 
spend the majority of their time in the Hard State (HS), 
which is associated with a low accretion rate, a hard X- 
ray power law component with photon index 1.4s; F s;2.1, 
and steady radio emission with an optically thick, flat-to- 
inverted spectrum. Sources that persist in the HS for several 
weeks generally show clear correlations between the radio 



and X-ray luminosities {Lr oc Lx) { Corbel 



et al. 2003 Gallo, Fender, fc Pooley 20031. More recently. 



2000 



Corbel 



evidence for a similar correlation between the near-infrared 



(NIR) and X-ray bands has also been found (Russell et al, 
2006{[Con"at et al.|200"9l ). 

Similar to the other BHBs, GRS 1915-1-105 also dis- 
plays periods of relatively hard, steady X-ray emission, but 
in contrast to the HS, only for about half the observation 
time ( |Trudolyubov|2001[ ). [Belloni et al.| ( |2000[ ) identify these 
intervals of decreased variability with long C state episodes. 
The subset of state C/ class x observations that are radio- 
bright are referred to as the plateau state (ibid.), but else- 
where have variously been been referred to as the radio-loud. 



radio-plateau or radio loud low/hard X-ray state (Muno 



et al. 20011, the type II hard steady state (Trudolyubov 



2001[ ) and xrl ( [Naik fc Rao|[2000 1 



The plateau state was first described by "Foster ct al.l 



1 1996 1 and its description was later refined in Fender et al.| 
(19991. The state can assert itself in a period as short as 



a day and is characterised by a decrease in X-ray flux den- 
sity and an increase in radio flux density to a steady value, 
typically ~ 10 — 100 mjy. Ample evidence supports the fact 
that the plateau state X-ray and radio luminosities are in- 
deed correlated, with an increasing delay from X-ray to in- 
frared (IR) to radio emission (e.g. Klein- Wolt et aL][2002 |. 
As in the HS the plateau radio spectrum is optically thick, 
and AU-scale self-absorbed compact jets have been spatially 
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Figure 1. HID comparing our observations (filled circles) to the 
plateau states from [Belloni et al.| | |2000| (data courtesy of T. 
Belloni; open stars are xl, closed stars are x3). The soft, higher- 
luminosity red dot is ObsID 90105-07-02-00, the lower is ObsID 
40403-01-09-00. The plot shows the hardness and fiux, normalised 
to Crab units, calculated by taking the ratio and the sum of the 
RXTE PCU2 rates in two bands: channel 0-35 (~ 2 - 15 keV) 
and 36-49 (~ 16 — 21 keV). Each band is divided by the Crab rate 
in the same band, on the same day, simulating a photon index of 
2.1 and utilising a normalisation of 10 photons cm~'^ s"'^ keV~^ 
at 1 keV. 



resolved using Very Long Baseline Interferometry (VLBI) 
( [Dhawan, Mirabel, fc Rodriguez[[2000| [Fuchs et al.[[2003b[ ). 
Another marker of the plateau state comes from timing anal- 
yses, where 1-10 Hz quasi-periodic oscillations (QPOs) are 
present ( |Rao et al.[2000| ). The exact frequency of these QPOs 
appears to be inversely correlated with the radio and soft X- 
ray flux ( [Ran fc Greiner|2003| [Rodriguez et al.[[2002 l. 

In line with these arguments, it is tempting to iden- 
tify the plateau state as GRS 1915+105's analog of the HS. 
However, although the plateau state and the HS share many 
similarities, it does display some distinct properties that can- 
not be ignored. For instance, while the BHBs in the HS 
usually have a luminosity of < 10% //Edd, the average lu- 
minosity observed in the plateau state is ~ I/Edd- Moreover, 
the plateau X-ray photon index is never very hard, with F ~ 
1.8 - 2.5 ([Fender fc Belloni|2004[ ). Finally, [Reig, Belloni, fc 
[van der Klis[ ( [2003[ ) argue that the canonical HS is never 
seen in GRS 1915-1-105 because a power-law tail (with no 
exponential cut off) is always present in the plateau state 
hard X-rays. Although the origin of such tails in the softer 
states of canonical BHBs is still under debate, it is generally 
not associated with the HS. 

In this work, we seek to make a more quantitative com- 
parison between the plateau state found in GRS 1915+105 
and the HS found in more typical black holes, in the con- 
text of an outflow-dominated model that has successfully 
described broadband data from several BHBs. In Section 
[2|we describe the multi- wavelength data and the reduction 
methods. We discuss what physical parameters we use for 
modelling and why in Section [S] and the model itself, to- 
gether with the results in Section [4] We put our work into 
context with the findings of others in Section[5] before draw- 
ing our final conclusions in Section [6] 
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Table 1. Observations included in dataset 1 and 2, with RXTE ObsID 40403-01-09-00 and 90105-07-02-00 respectively. 
The H and K bands, obtained using the Cooled Grating Spectrometer at UKIRT, are at 1.455-2.094 and 1.906-2.547 ^m, 
respectively (J band measurements were also done, but these are unusable due to high interstellar extinction). The K band 
observation in dataset 2 was done (more than half a day before ObsID 90105-07-02-00) with the ANDICAM at CTIO. The 
GBI measurements are 2.25 & 8.3 GHz and the flux densities used are 30 sec vector scan averages. The Ryle telescope operates 
at 15 GHz and the flux density used is the average of ^13 ks worth of five minute integrations. 



Dataset 1 1999 July 8 / MJD 51367 



Dataset 2 2005 April 13 / MJD 53473 



Band 



Instrument 



MJD start 



UTC 



Duration 



Instrument MJD start 



Time 



Duration 



X-ray 
NIR 

Radio 



RXTE 




51367.2737 


UKIRT 


{H) 


51367.4821 




(K) 


51367.8216 


GBI 




51367.371 



06:34:13 14435 sec 
11:34:10 80 min 
07:56:59 64 min 
09:02:53 30 sec 



RXTE 53473.054 01:09:05 6769 sec 

CTIO [K) 53472.3216 07:43:06 25 min 



Ryle 



53473.098 02:21:07 219 min 
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Figure 2. GRS 19154-105 X-ray and radio lightcurvcs for the 
1999 data set. Also the radio spectral index, as measured by the 
Green Bank Interferometer is shown. The red dashed line indi- 
cates the time of the RXTE observation. 



•I CD 140 

1^120 
^^100 



CO 



CO 






CD 
•D 

3 



CO 



CD 

D 



80 
60 

40 

13 

12.8 

12.6 

12.4 

12.2 

12 

11.8 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 





^^-w^^* 



ASM Lightcurve h 



*3*^ 



± 



i=— \ 



K band CTIO 



± 



15 GHz Ryle 



%:* 



.+* 



^4 






^+a 



A* 



i:*^ 



•tiH-+ 



53420 53440 53460 
MJD 



53480 



Figure 3. GRS 1915-fl05 X-ray and radio lightcurves for the 
2005 data set. The radio lightcurve was measured by the Ryle 
Telescope jPooley fc Fender|1997^ . The red dashed line indicates 
the time of the RXTE observation. 



2 OBSERVATIONS AND DATA REDUCTION 

Based on the lightcurves and its position in the Hardness 
Intensity Diagram (HID) diagram (see Figures [l] |2] & |3| , 
GRS 1915+105 was in a class x state on July 8th 1999 
(dataset 1) and April 13th 2005 (dataset 2). The presence 
of a 2-5 Hz QPO on both occasions (see Figure |4]) corrobo- 
rates this fact. During the 1999 episode, quasi-simultaneous 
radio and near-infrared observations were performed using 
the United Kingdom Infra-Red Telescope (UKIRT) and the 
Green Bank Interferometer (GBI). The 2005 dataset also 
holds observations of the Cerro Tololo Inter-American Ob- 
servatory (CTIO) and Ryle telescope. Details on the obser- 
vations and reduction methods are found in the individual 



Sections and Table [T] The 2005 RXTE data is also exten - 
sively analysed and discussed in Rodriguez et al. ( 2008a|b I 



2.1 X-ray: RXTE data reduction 

We use data from two instruments on board the RXTE: 
the Proportional Counter Array (PCA; [Jahoda et al.|[2006[ ) 
and the High Energy X-ray Timing Experiment (HEXTE; 
|Rothschild et alT]|1998 l. The X-ray data have been reduced 
using HEASDFT, version 6.3.1., applying standard extraction 
criteria: a pointing offset of < 0.01° from the nominal source 
position, and a source elevation of > 10° . The exclusion time 
for the South Atlantic Anomaly (SAA), however, was only 
10 min, as GRS 1915-fl05 is a very bright source. For the 



4 P.C.N, van Oers et al. 
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100 



Figure 4. Normalised (following Leah y et al.||1983l l power den- 
sity spectra showing the type-C QPOs in the X-rays for both 
datasets, indicating the observations are plateau state. (Courtesy 
of P. Soleri) 



same reason we disregarded the top layer and imposed a 
maximum "electron ratio" of 2, to take into account con- 
tamination by source photons. All HEXTE data products 
were dead-time corrected. Correction for the PCA dead-time 
was also carried out; PCU 4 was inactive for half of the 
1999 observation, while during the other half, both PCU 1 
and 4 were inactive. The 2005 observation was done with 
PCU and 2 only. Due to uncertainties in the first four 
PCA channels, we only include PCA data above 3keV. The 
PCA calibration falls off rapidly above ~ 25 keV and since 
HEXTE provides reliable data for energies > 20 keV we ig- 
nored the PCA data above 22keV. For the PCA, standard2f 
mode data were used. The PCA background was modelled 
using the pca_bkgd_cmbrightvle_eMv20051128 model. Only 
HEXTE data above 20 keV were used due to the uncertainty 
of the response matrix below these energies. At high ener- 
gies, the spectrum was ignored above 200 keV. We rebinned 
the X-ray data to a minimum S/N of 4.5. 
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Figure 5. NIR spectrum of the 1999 data set 

librarjrl resampled to the CGS4 resolution and flux-scaled 
to match the HD 176232 magnitudes in H and Ks. The un- 
certainty of the flux calibration obtained this way is about 
5%. 

The spectra were then dereddened using the extinction 
laws given in ( Cardelli, Clayton, fc Mathis|1989p and ( Chiar 



& Tielens 20061, assuming a hydrogen c olumn density of 
A^H=3•5xlF^m-^ i.e. Av=19.6-h/-1.7 ([Chapuis fc Cor- ' 



|bel||2004l The errors resulting from using the two different 
methods were then added in quadrature. The overall prop- 
agated uncertainties, combining flux calibration and dered- 
dening, represent about 20% of the flux. The final spectrum 
is shown in Figure [5] 

The UKIRT spectrum included in dataset 1 was also re- 
binned, because the spectral resolution is well beyond what 
we need for continuum fitting (see Section B|. 

The 12.40 mag K band point in dataset 2 was taken 
using the CTIO 1.3m telescope using the ANDICAM (A 
Novel Double-Imaging CAMera) detector (Neil, Bailyn, & 



Cobb 20071. It was dereddened using the same method as 



the UKIRT data. However we chose the error bars to account 
for the fact that GRS 1915-1-105 may have been approaching 
a state transition: The source magnitude magnitude varies 
from 11.8 to 12.72 in about 9 days (see Figure jsl, suggest- 
ing a dereddened flux density of 67.85±27.25 mjy. On MJD 
53473 we also see the flux in radio increasing (see below). 



2.2 NIR: UKIRT data reduction and CTIO 

On 1999 July 08, medium-resolution spectroscopic obser- 
vations of GRS 1915-1-105 were performed with the Cooled 
Grating Spectrometer (CGS4) mounted on the United King- 
dom InfraRed Telescope (UKIRT, P.I. Harlaftis). Twenty- 
four spectra in the Bl and eighteen in the B2 filters were 
acquired, each time at two nodding positions, giving a total 
exposure time of 1260 s in Bl and 1080 s in B2. Moreover, 
the telluric standard HD 176232, a FOp main sequence star, 
was observed in the same conditions. 

All the spectra were reduced with IRAF through bad 
pixel correction, fiat fielding, and sky subtraction. The tar- 
get spectra were then extracted and wavelength calibrated 
by extracting, in the same condition, krypton and argon 
spectra in Bl and B2, respectively. They were combined and 
telluric feature-corrected through the division by the stan- 
dard star spectrum. We finally multiplied the corrected spec- 
trum by a FOV star one from the ISAAC synthetic spectra 



2.3 Radio: GBI and Ryle Telescope 

Dataset 1 includes two GBI data points and dataset 2 in- 
cludes a data point from Ryle Telescope. Lightcurves from 
these telescopes are shown in Figures [2] and [31 respectively. 
For the GBI data, a fiux density calibration procedure 
similar to that reported in Waltman et al (1994) has been 
employed here. The fiux densities of 0237-233, 1245-197, and 
1328-1-254 were determined using observations of 1328-1-307 
(3C 286). The fiux density of 3C 286 was based on the scale 
of Baars et al. ( 1977 1, and the assumed values were 11.85 Jy 
at 2.25 GHz and 5.27 Jy at 8.3 GH^ On MJD 51367.371 
the GBI fiux densities were 25±4 mJy at 2.25 GHz and 
30±6 mJy at 8.3 GHz. Similar behaviour was shown during 
a plateau state of April 2000 (lUeda et al.||2002l). 



^ http : //www . eso . org/sci/facllities/par anal/instrument s/lsaac/tools/lib/i 
^ http : //www . gb . nrao . edu/f gdocs/gbl/plgbi/README 
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The Ryle Telescope operates at 15 GHz and it observed 
GRS 1915+105 (simultaneously with ObsID 90105-07-02- 
00) from MJ D 53473.098 to MJD 53473 .417. The data used 
here are from Rodriguez et al. (2008a I and were reduced 



following Pooley & Fender (19971; observations of Stokes 



I-l-Q were interwoven with those of a nearby phase calibrator 
(B 1920-1- 154). The flux-density scale was set by reference to 
3C48 and 3C286, and should be consistent with that defined 



by |Baars et al.| ( |1977 |. 

On MJD 53473 we see the average radio flux density 



increasing from 44.9±3.0 to 70.4±4.3 mjy (Rodriguez et al 



2008b I, in addition to the IR lightcurve also suggesting that 



in spite of the HID classification, GRS 1915+105 may not 
have strictly been in a plateau state. For this dataset we 
therefore use the average from MJD 53473.098 to 53473.25 
of 44.9+3 mJy, closest to the RXTE and CTIO observations. 



3 CONSTRAINTS ON INPUT PHYSICAL 
PARAMETERS 

In order to meaningfully compare broadband spectra, we 
will fit our datasets using a model designed for simultane- 
ous radio through X-ray datasets. Historically the modelling 
of BHBs has focused on accretion flow models of only the 
X-ray emission, of which Comptonising corona models have 
been particularly successful. However in recent years, the 
evidence is mounting that the bipolar jet outflows found in 
the HS contribute significantly across the broadband contin- 
uum. For bright transients, e.g. jRussell et al.| ( |2006[ |2010[ ) 
estimate up to 90% of the NIR and up to 100% in the X-rays 
could be dominated by jet production. Other observations of 
radio/IR/X-ray correlations ( [Corbel et aLlpOOO) [2003} [Co] 
|riat et al.|2009 1 can be interpreted by either synchrotron or 
synchrotron-self Compton (SSC) processes relating to the 
jets. The outflow-dominated model of|Markoff, Nowak, fc| 



Wilms ( 2005 1 (hereafter: MNW05) is still the only model 



that can fit the X-ray spectrum with the same precision as 
corona-only models (ibid; Nowak et al., in prep.) while also 
fitting the radio through IR bands from the same physical 
picture. The MNW05 model has already been applied suc- 
cessfully to many Galactic sources in the "canonical" HS: the 
original paper features fits to Cyg X-1 and GX 339-4, while 
a different data set of the latter, along with observations 
of XTE Jl 118+480 are fit in jMaitra eFal] ( [2009| . Further 
papers have explored applications to simultaneous broad- 



band data sets from GRO J 1655-40 and A0620-00 ( [MigUari 
[et al.|2007|[Gallo et al.|2007| l, with other sources in progress. 
The results of all these applications has been the discovery, 
perhaps not surprisingly, that the free parameters fit into 
similar ranges for all stellar mass sources. 

Furthermore, the discovery of a Fundamental Plane of 



black hole accretion (Merloni, Heinz, & di Matteo 2003; Fal- 



cke, Kording, & Markoff 20041 supports mass-scaling ac- 



cretion physics from stellar to supermassive BHs, and thus 
would argue that the MNW05 should also apply to weakly 



accreting AGN. Confirming this, MarkofF et al. (20081 suc- 



cessfully fit several spectral energy distributions (SEDs) 
from the supermassive BH M81* with parameters in the 
ranges found from BHBs. 

There exists, therefore, a solid framework for the mod- 
elling of BHBs in the canonical HS, against which to test 



fits of data from the GRS 1915+105 plateau state. In this 
paper we do not consider whether MNW05 is correct but 
rather use it as a benchmark to test how far GRS 1915+105 
fits the phenomenology of the hard state. 

Whenever possible, we fixed values for model physical 
parameters in accordance with those found in other publica- 
tions during the fitting process. In the following subsections, 
we first discuss how we obtained these values (presented in 
Table [2|, and then briefly summarise the MNW05 model 
and our fitting methods before presenting our results. 



3.1 Distance and Hydrogen Column Density 

The distance to GRS 1915+105 is still a matter of some de- 
bate. A first estimate was derived in [Mirabel fc Rodriguez[ 
( jl994 |. From a core ejection they find a maximum distance 
of 13.7 kpc, under the assumption that the ejection was in- 



trinsically symmetric. Rodriguez et al. ( 1995 1 attempted to 



derive a more accurate distance estimate by determining the 
kinematic distance from 21 cm absorption spectra of atomic 
hydrogen along the line of sight to GRS 1915+105 during 
a radio outburst, and find that it could be as far away as 
12.5+1.5 kpc. Later measurements of the ^^CO(J = 1 — 0) 
spectrum by Chaty et al. ( 1996 1 are consistent with this 
distance. They also constrain the visual extinction to be 
Av = 26.5 ± 1.7 mag, corresponding to a total hydrogen 
column density of A^h = 4.7 + 0.2 x 10^^ cm"^ 

If GRS 1915+105 resides at 12.5 kpc, it must be accret- 
ing above the Eddington limit (Lx-ray ~ 2.9 x 10'^^ erg s~^; 
McClintock & Rcmillard 2003 giving Lx-ray /iEdd ~ 1.5, 
using a mass of 14 M©; see next Section). A slightly lower 
maximum distance of 11.2+0.8 kpc, derived from proper mo- 



tion studies of radio ejecta during four major events (Fender 
[et al.|1999| ), would still give a ratio Lx-ray /^Edd ~ 1.3. 

Ho wever, this distance could be too large. [Chapuis fc[ 



Corbel (20041 measure CO(J =1-0) velocity spectra 



of clouds along the line of sight, as well as of two nearby 
HII regions, and re-evaluate the hydrogen column density to 
Nn = 3.5 ± 0.3 X 10^^ cm~^. This column gives a visual ex- 
tinction of A„ = 19.6+ 1.7 mag, reducing the number inter- 
vening molecular clouds, and arguing for a smaller distance 
to GRS 1915+105 of 9.0+3.0 kpc (at 9 kpc the accretion rate 
estimate drops below Eddington), based on evidence that 
GRS 1915+105 lies behind the molecular hydrogen cloud 
G 45.45 +0.06. This cloud is at ^7 kpc ( [Feldt et al.[[1998[ ), 
however [Chapuis fc Corbel] ( |2004| suggest that this cloud 
belongs to a huge complex, located at a distance of 6 kpc, 
providing the above lower limit to the distance (the upper 
limit of 12 kpc is from [Fender et al.[|1999[). However, the 



high column density required (Mirabel & Rodriguez 1994 



Rodriguez et al. 1995 1 suggests that another molecular cloud 
should be between us, therefore jFender et al.[ ( [T999| ) adopted 
the conservative distance of llkpc. 

Considering the assumptions made in [Chapuis fc Cor-[ 
bel[ (|2004| and the resulting size of the errorbars, we have 



decided to follow Fender et al. ( 1999 1 in adopting a more 



conservative distance of 11 kpc for this work. This distance 
is consistent with Nh ~ 3.5 + 0.3 x lO^'^ cm~'^, which we be- 
lieve is the correct value pertaining to the IR band, as well as 
the larger hydrogen column density of 4.7x10^^^ cm~'^ from 
[Chaty et al.| ( [1996| ), that we use for our X-ray fits, motivated 
by the works of e.g. | Klein- Wolt et al^ ( [2002[ ) (see Section[5|. 
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Table 2. Fixed physical parameters used in the fitting process, 
obtained from the literature. We omit the error bars, as they are 
not used in the fits. 



parameter 



value 



units Reference 



column density 

mass 

inclination 

distance 

donor temperature 



4.7 

14 and 23 

66 and 50 

11 and 6 

4455 



102 



cm ^ 


a 


M0 


b 


o 


c 


kpc 


c, d 


K 


e 



References; a: |Chaty et al.| ||1996[l, b: [Greiner et al.| | |2001a[l, c: 
Fender et ah] ||1999|l, d: |Chapuis fc Corbel| l |2004| |, e: |Alonso,| 
Arribas, fc Martine2-Roger|||1999[l 



Because a distance of 11 kpc implies super-Eddington 
luminosities for GRS 1915+105, we also fit the data using 
the minimum distance of 6 kpc, in order to explore the im- 
portance of distance on the modelling conclusions. 



3.2 GRS 1915+105 Black Hole Mass 

In order to constrain the mass of the compact object in a 
binary via the mass function, the orbital period, inclination 
and the mass of the companion (or the mass ratio) must all 



be known. Greiner, Cuby, & McCaughrean (2001a) found 
an orbital period Porb of 33.5 days and a velocity amplitude 
_ft'=140±15kms"\ giving a mass function f{M) = 9.5 + 3.0 
M© . Assuming the binary plane is the same as that of the 
accretion disc, the orbital inclination can be estimated from 
the orientation of the jets. As no constant precession has 
been observed, the jets can be assumed perpendicular to 
the accretion disc and orbital plane. The exact inclination 
is however still open to debate, as it is determined from 
the brightness and velocities of both the approaching and 
receding blobs. At a distance of 11 kpc the inclination is 
i = 66°±2° ( [Fender et al.|1999| |. [Ha"rlaftis fc Greiner] ( [2004 1 
were able to deduce a mass ratio Md/Mx = 0.053 + 0.033 
from the rotational broadening of photospheric absorption 
lines of the donor star. Together these values give a mass 



for GRS 1915+105 of 14.0±4.4Mq (Greiner, Cuby, fc Mc- 



Caughrean 2001a I. We use these as our fiducial values for 



the distance and mass. For the 6kpc distance we calculate 



a smaller inclination of i = 50° using Fender et al. (19991, 



and, using Harlaftis & Greiner (20041 a mass of 23M( 



3.3 Properties of the Donor Star 

A rough identification of the companion was performed by 
Greiner et al. (2001b I, by analysing absorption line fea- 



tures in the NIR. They conclude that the companion is a 
class III K-M giant and compare it successfully to a K2 III 
star. If the star is indeed a K2 III giant, its temperature 
should be T j = 4455+190 K ( [Alonso, Arribas, fc Martinez-[ 
Roger 19991, and accretion proceeds via Roche- lobe over- 
flow ( Greiner et al.|2001b l. We use this temperature for the 
simplistic blackbody spectral component representing the 
companion in our spectral fits. 



4 MODELLING AND RESULTS 



4.1 Description of model 



For all spectral fits we use the program ISIS (Houck & Deni- 
|cola[|2000[) com piled with XSPEC version 12.3.1x libraries 
( Arnaud|1996 l. The model discussed below is forward-folded 
through the X-ray detector response matrices, but applied 
directly to the radio through NIR data. To account for the 
additional uncertainties in the PCA response matrix, a 0.5% 
systematic error has been added in quadrature to all PCA 
data. As the relative calibration of the PCA and HXTE in- 
struments is not certain, the normalisation factor for the 
PCA data is set to unity and tied to the radio and NIR 
normalisations during the fits, while the HEXTE data nor- 
malisation is left free to vary. 

As the focus of our paper is on modelling the non- 
thermal spectrum, and the exact stellar model needed for 
GRS 1915+105 is uncertain, we use a simple blackbody to 
model the NIR, with the temperature fixed as discussed in 
the previous Section. This component serves to account for 
the excess IR/optical flux level due to the star, so that the 
overall model normalisation is correct. 

In addition to the jet and the star we have an outer 
accretion flow in the form of a "standard" geometrically 



thin accretion disc (Shakura fc Sunyaev 19731. This flow 



is parametrised by the radius of the inner and outer disc 
edges (rin,rout) and the temperature at the inner edge (Tin). 
In the Schwarzschild geometry, the innermost stable circular 
orbit has a radius of 6 r^, but this radius can be reduced to 
1 Xg in the Kerr metric, when the BH is maximally spinning. 
The disc component serves to both fit spectral signatures of 
thermal emission from the accretion flow in the X-ray band, 
as well as contributing seed photons for inverse Compton 
scattering, but is only a weak element in comparison to the 
entire spectrum. 

We now summarise the physical parameters in our jet 
(see Table [sl . One of the most important of these parame- 
ters is the jet luminosity {Nj). This factor scales with the 
accretion power at the inner edge of the accretion disc and 
represents the power going into the jet and acts as a nor- 
malisation factor. The power is equally divided to supply 
the internal and kinetic pressures. We assume that the ki- 
netic energy is carried by cold protons, while the leptons do 
the radiating. The energy involved in the internal pressure 
goes into the particle and magnetic energy densities, with a 
ratio determined by fc. fc = 1 equals equipartition and higher 
values indicate magnetic dominance. The initial velocity at 
the base of the jet, or "nozzle" is the proper sound speed of 
an electron/proton plasma: /37c ~ 0.4c. The radius of the 
jet-base is also a free parameter, ro (expressed in in units 
of gravitational radius Vg = GM/c?). The particles start of 
in a quasi-thermal (relativistic) Maxwellian distribution, the 
peak energy of which is determined by the electron temper- 
ature (Te). Beyond the nozzle the jet is allowed to expand 
freely, or adiabatically, causing a longitudinal pressure gradi- 
ent. This pressure gradient leads to a moderate acceleration 
of the jet along the direction of flow; the resulting velocity 
profile is calculated from the Euler equations and is roughly 
logarithmic, evening out at bulk speeds with Lorentz factors 
of r ~ 2 - 3 (see e.g. [Falcke, Markoff, fc Bow[[2009[ ). In a 
segment of the jet located at some variable distance (zacc; 
also expressed in rg) from the base, we assume a significant 
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fraction (75 %) of the leptons is accelerated into a power- 
law. The slope of the power-law particle distribution p is also 
not pre-determined. In each jet segment after the accelera- 
tion front, the general shape of the distribution is assumed 
to remain the same, while the total lepton density decreases 
according to the adiabatic expansion. This is achieved as- 
suming a continuous injection of "fresh" energetic power-law 
distributed electrons {N{E)dE oc E~'^dE) in each segment 
after the acceleration region. The jet is stopped at a distance 
Zmax from the base. 

All fits are done using the following components: (1) 
The MNW05 steady-state outflow-dominated model that in- 
cludes a multi-colour blackbody accretion disc and a single 
blackbody for the companion star; (2) an additive Gaussian 
line profile, with a line energy left free to vary between 6 
and 7 keV, and line width a free to vary between and 
2 keV; Models (l)-|-(2) are either convolved with Compton 
reflection from a neutral medium (reflect; Magdziarz & 



|Zdziarski|1995 1 or multiplied with the smeared edge model 
(smedge; |Ebisawa||1991 using the "standard" index for the 
photoelectric cross-section of ~ —2.67;), that accounts for 
relativistic smearing of the iron line, and multiplied with a 
photo-electric absorption model (phabs) to account for the 
interstellar medium. As the strength of an absorption feature 
(or edge) is related to the strength of the according emission 
line, the line width is in principle also related to the absorp- 
tion edge width. Thus when using the smedge model, we tie 
the width of the edge to the width of the Gaussicoi. For the 
reflect model we assume the viewing angle to correspond 
to the jet inclination. Furthermore, although the reflect 
model already includes an absorption edge, it is a sharp un- 
smeared edge, therefore we also tried fitting the data with 
both the smedge and reflect models. 

4.2 Spectral fitting and results 

The results of all fits are presented in Figures lOHll] with 
respective parameters listed in Table [3] Below we discuss 
the individual fits to both datasets in more detail. 

Initial fits revealed that, in the case of GRS 1915-f 105, 
using smedge is statistically preferred to using reflect, to 
model the refiection features above ~ 10 keV. Therefore, 
we proceed to fit the data in greater detail with the smedge 
model, but include a fit that employs reflect as this model 
has been used generally to model the disc reflection in pre- 
vious works (see Section [s]). 

The spectral index of the optically thick radio-NIR syn- 
chrotron emission is determined in part by the internal jet 
plasma parameters such as the electron temperature, but it 
is most sensitive to the Doppler beaming factor (calculated 
from the inclination of the jets). The closer the jet axis aligns 
with the line of sight, the less inverse the observed radio 
spectrum. Fixing the inclination according to observation 
and choosing a jet length for all fits of 10^^ cm, sufficient 
so that the slope through the radio data points is continu- 
ous, reveals a key difference between the two GRS 1915+105 
plateau state datasets and canonical BHBs. In order to avoid 
too much excess in the NIR over the companion star BB, 
and to provide the best NIR fit, the initial particle accel- 
eration region in the jets Zacc needs to be at a distance of 
at least ~ 10* Tg. Therefore the post-acceleration region in 
the jets does not dominate the X-ray emission below ~ 10 



keV, in contrast to canonical BHBs in the hard state where 
acceleration is found to begin on the order of 10s of r-g . 

In general we obtain values for electron energy index 
p ~ 2.2 that result in a significant synchrotron contribu- 
tion to the X-rays below ~ 20 — 50 keV, with a ratio of 
synchrotron to inverse Compton flux of ~ 0.1. The steeply 
declining X-ray flux above ~ 20 keV is best fit by a domi- 
nant inverse Compton contribution from the base of the jets. 
With a much harder value of p we could conceivably fit more 
of the X-ray emission via synchrotron emission, but only if 
the exponential decay shape plays a significant role (see Fig- 
for such a fit to dataset 2). For very soft values of 
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the spectral index, synchrotron emission will not contribute 
significantly, but fits without any synchrotron contribution 
to the X-ray at all are not statistically favoured. 

The exact spectral shape and normalisation in the post- 
acceleration synchrotron component (fitting the radio data) 
are very dependent on the value of p. When calculating a lo- 
cal statistical minimum, the fitting routine naturally favours 
the plethora of X-ray data points over the few radio data 
points, in finding the confidence limits for p, explaining non- 
negligible systemic residuals in the radio. 

For clarity we note that below we will use SMALL CAPS 
to refer to fits and will continue to use the typewriter font 
to refer to models. 



4.2.1 Dataset 1 

The best fit models for dataset 1 are shown in Figures [6] - 
[9] (corresponding to Table [s] smedge, reflect, smedge2 
and 6kpc column respectively; see below for explanation of 
these fit names). The fits are in general statistically good, 
with a large contribution to the residuals coming from the 
inadequacy of a simple black body to fit the NIR data (see 
Figure |5|. Completely removing the NIR data from the best 
fit SMEDGE models an improvement in the AXrcd of ~ 0.08, 



to x?cd 



0.69. 



Our best fit smedge, employing the smedge model also 
shows the most extreme behaviour in terms of the distribu- 
tion of the energy budget (evident from the value for the 
magnetic dominance k) and the acceleration front distance 
Zacc- In addition this fit is likely to suffer from pair produc- 
tion (see TableWland Sectionlsl. Hence we explore several fits 
with less extreme physical properties, including two more fits 
employing a smedge model, (smedge2 and 6kpc) and one 
fit using the reflect model (reflect). The smedge and 
SMEDGe2 fits are different local minima that primarily dif- 
fer from each other in the distribution of the energy budget, 
most importantly the latter fit has a reduced value for the 
magnetic dominance k and particle distribution index p. For 
6kpc we reduced the distance to GRS 1915-1-105 to 6 kpc 
(see below). These fits clearly show over an order of mag- 
nitude decrease in z^cc and a factor 2 — 25 reduction in k. 
Moreover, pair production will be less important. We note 
here that fixing Zacc to 2 x 10* rg in the smedge fit results in 
only a minor increase in x?cd of 0.05, while producing virtu- 
ally the same model parameters. This suggests that in this 
case our model does not converge and shows artefacts from 
the fitting routine that tries to find the absolute local mini- 
mum (within the tolerance limit), although there is no real 
physical basis for this insignificant improvement. Hence for 
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Table 3. Parameter ranges found in canonical black holes GX 339-4 l|Markoff et al.|2005[ Maitra, Markoff fc Bro cksopp 2009), XTE 
J1118+480 (Maitr a, Markoff & Brocksopp 2009), Cyg X-1 ( Markoff et al.|2005| l, GRO J1655-40 ( [Migliari et"alf20 07) and A0620-00 
jGallo et al.|2007| , during the HS and best-fit parameters for the GRS 1915+105 plateau state. The error bars have been resolved 
at 90 percent confidence level. We failed to resolve the error bars for parameters listed in italics. A^^ is the jet normalisation, ro 
the nozzle radius, Te the temperature of the leptons as they enter the jet, p the spectral index of the radiating particles, k the 
ratio between magnetic and electron energy densities, Zacc location of the particle acceleration region in the jet, /scat 0.36/ ratio 
between scattering mean free path and gyroradius, L^isc a^nd T^^^^ the luminosity and temperature of the accretion disc, Ahxt is 
the relative normalisation between HEXTE and PGA, Aw^^,, inline and a are the normalisation, the energy and width of the iron 
line feature, and Sedge and Tmax are the energy and the normalisation of the smedge model iron line edge (the iron edge width is 
the same as that of the Gaussian) and Q,/2-k is the reflection fraction from the reflect model. 



variable 



units 



range found 

in canonical 

BHBs 



GRS 1915+105 



MJD 51367 



SMEDGE2 



6kpc 



MJD 53473 



''0 
P 

k 



.^disc 
-^ disc 
^HXT 



lO-^LEdd 

GM/c^ 
10^ K 



lO^GM/c^ 

10-4 

GM/c^ 

10-1 LEdd 

keV 



0.0034 - 0.71 

3.5 - 20.2 
20 - 52.3 
2.1 - 2.9 

1.1 - 7 
0.007 - 0.4 

1.6 - 299 
0.1-486 

0.007 - 0.99 
0.06 - 1.53 



4.78"' 



;i+0.01 
-0.00 
20 39+001 
q 21+0.00 
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^•■-iu_o.o6 

692^1 
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6.45^ 



-+0.00 
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9.9r 
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0.00 
+0.01 
+0.01 
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2.22 
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+0.00 

0.00 
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i-^^-o.oo 

+0.0 
0.1 

r, Q7+0.02 

4 67+0-04 
^•°'-0.06 

1.51 

+0.01 
0.00 



+0.01 

-0.01 

4 q4-t-0.05 

^--'^-o.oi 
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0.04 

+2 
3 
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^+0.00 
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707 



+68 
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+0.03 

0.04 
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113/147 
(=0.77) 
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(=1.02) 


166/147 
(=1.13) 


131/147 
(=0.89) 


59/81 
(=0.73) 


59/81 
(=0.73) 



MNW05 Canonical Parameter Ranges and GRS 1915-|-105 Best-Fit Parameters 

t Value is pushing lower or upper boundary and is therefore not well constrained, 
b Derived from model values for r^^ and T^iac- 



the SMEDGE fit we will ignore the value for Za, 
the other parameter values in our analysis. 



but include 



The steeper optically thin synchrotron spectrum of 
SMEDGE2 could be reconciled with the steeper p found in 
canonical BHBs by assuming that both share such a hard 
injected power-law, and considering the effect of cooling on 
the spectrum. We do not consider cooling explicitly here, 
but are investigating such issues in a separate work. Because 
of the diminishing magnetic field strength along the jets, 
the cooling (dominated by synchrotron losses) is negligible. 
The break energy where the particle distribution steepens 



in "particle space" shifts according to E oc B our (Kar- 



dashev|1962 l. Thus the further along the jet we go, the less 
important cooling becomes, and thus the largest contribu- 
tion to the higher frequencies from synchrotron radiation is 
from the first acceleration zone. Because this zone occurs al- 
most two orders of magnitude further out in GRS 1915+105 
compared to BHBs, this effect could explain why we are see- 
ing the uncooled injected spectrum rather than the cooled 
spectrum, expected to be steeper by 0.5. 



is clear that the true distance to GRS 19154-105 could be 
much smaller than the 11 kpc adopted for most of our fits. In 
an attempt to explore the effect of distance, we performed 
one fit on dataset 1 with the distance to GRS 1915-1-105 
reduced to 6 kpc (see Figure |9] and Table [s] 6kpc column). 
The modifications result in a similar x?ed but some values 
are closer to what we have come to expect for the HS in 
other black holes. Firstly we find Zacc to be the smallest of 
all fits. This can be understood from the smaller inclination 
of 50° corresponding to the smaller distance, resulting in 
a less inverse radio spectrum and hence a lower flux level 
at the IR break. In addition the jet luminosity Nj and the 
partition parameter k are closer to values found for other 
BHBs. With Nj ~ 0.22LBdd and A: ~ 90 they are only a 
fraction of the values found for most 11 kpc fits. Also the 
smaller nozzle radius ro ~ 5 is more typical of the usual 
HS value. One possible interpretation of these results is that 
GRS 1915+105 might in fact be closer than the conservative 
value usually taken in the more recent literature. 



From the arguments outlined in the previous Section it 
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4.2.2 Dataset 2 

The best fit models for dataset 2 are sliown in Figures p^ and 
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(corresponding to Table [s] column, eltemp and SYNCH 
respectively), eltemp offers the best comparison - in terms 
of what component fits what feature - to dataset 1, and owes 
its name to its low electron temperature Te, while SYNCH is 
dominated by the post acceleration synchrotron component. 

Fitting dataset 2 with similar parameter values (relating 
to the energy budget and assumed geometry) as dataset 1 
yields very poor statistics (best xicd ~ 30). Trying to model 
the steep X-ray spectrum, mainly employing the exponential 
decay of the multicolour disc fails, because a higher disc 
contribution offers more (soft) seed photons. These photons 
are up-scattered and create a Compton tail in the hardest 
part of the spectrum where they overestimate the observed 
flux. 

We therefore searched the parameter space for any other 
solution that reduces the flux of Comptonised high-energy 
photons by reducing the electron temperature Te. A lower 
temperature electron distribution would on average not up- 
scatter the disc seed photons to equally high energy. We 
flnd that if we allow Te to evolve freely, we end up with a 
very good flt (x?ed ~ 0.7), however the final temperature is 
rather low. At ~ 4 x 10® K the bulk of the electrons would 
be sub-relativistic. The use of a relativistic Maxwell- Juttner 
distribution, as done in the MNW05 model, would no longer 
be fully justified, as the majority of the particles in the ther- 
mal distribution would have 7 = 1. However this approach 
effects a very peaked and steep pre-shock synchrotron com- 
ponent. This distribution of seed jet-base photons allows for 
a Comptonised photon distribution that is steep enough to 
fit the X-ray data up to 50 keV. Above 50 keV the Comp- 
tonised disc photons harden the model spectrum just enough 
to fit the entire range of HEXTE data. 

An alternate possibility for fitting dataset 2 - with an in- 
creased electron temperature - is found when we increase the 
ratio of synchrotron to inverse Compton emission. With syn- 
chrotron dominating below ~ 30 keV (see Figure ITT] and Ta- 
ble |3| SYNCH column) and cutoff at about 1 keV, the subse- 
quent exponential decay can approximate the X-ray features 
rather well, when combined with an accretion disc peaking 
at the same energy (Figure 111. In general we feel that, 



given the unknowns in the exact shape of the particle dis- 
tribution around the cutoff, it is undesirable to rely on this 
in general during the fits. However it is worth noting that 
it provides a very good description of the data (x?od ~ 0.7). 
Fitting with the synchrotron cutoff requires a very hard syn- 
chrotron spectrum (F ~ 1.2). Normally such a hard spectral 
index would be expected only in ultra-relativistic sources 
(e.g. |Heavens fc Drury|1988[ ). As for the SMEDGe2 fit, the rel- 
atively shallow spectrum could, however, again be (partly) 
reconciled with the canonical values for p, considering the 
effect of cooling on the spectrum. 



4.3 Stellar companion spectrum 

While the predominant temperature of a K2 star is 4455K, 
the dataset 1 fits show a systemic discrepancy that could be 
resolved increasing the blackbody temperature (and normal- 
isation). While the temperature of 4455 K is the predomi- 



Table 4. Relative importance of pair production processes at 
the base of the jet, for the statistically good fits in table [3] n is 
the lepton number density per cm~^, ripa and ripp are the pair 
annihilation rate and production rate per cm~^s~^, respectively. 
In order to compare the production rates to the lepton density, 
they have been multipUed by an estimate for the residence time, 
given by rO/c. If n 2> Max(npa,npp), and ripa > »ipp we should 
be safe from pair production. From these two considerations, the 
former is the most important. 



fit 


rO in T-g 


log(n) 


log(npa X rO/c) 


log(npp X rO/c) 


SMEDGE 


20 


15.3 




12.9 






15.0 


REFLECT 


9.3 


16.1 




14.3 






15.7 


SMEDGE2 


6.4 


16.6 




15.1 






15.8 


6kpc 


4.9 


16.0 




14.9 






15.1 


ELTEMP 


4.2 


17.1 




16.3 






13.4 


SYNCH 


3.6 


16.4 




14.5 






15.2 



nant value for a solitary class 2 star of this type, clearly this 
is not the case for the companion of GRS 1915+105. 

Although the physical properties of stars under these 
conditions are still largely undetermined, |Kaper| ( |2001[ ) re- 
solved discrepancies in the effective temperatures of donor 
stars in high mass X-ray binaries, of 10-25% higher than 
expected from their spectral classification. Increasing the 
temperature of our blackbody by 25% to 5500 K results in 
an insignificant improvement in fit by Ax?cd = 0.03. 

Letting the temperature and normalisation of the black- 
body in our best-fit result SMEDGE free to vary improves 
the Xrod to 81/165 = 0.5, but yields an extremely high 



blackbody temperature of 1.33 



+0.14 



X 10^ K. Such an ex- 



ceedingly high temperature indicates that other IR contri- 
butions than photospheric emission from the secondary are 
likely. |Chaty et al.| (1996) already concluded the same from 
the rapid variations and spectral shape they observed in this 
band. In particular irradiation of the companion star and/or 
the disc for a source as bright as GRS 1915+105 should 
be quite siginflcant, but the former is not accounted for in 
detail in our modelling, while the latter is not accounted 
for at all. Using an irradiated disc+jet model, |Maitra et al.| 
( 2009[ ) obtained disc temperatures for XTE J1118+480 in 
the order of a few tens of thousands Kelvin, in agreement 
with e.g. [Hynes et ah] ( |2006[ ). Our obtained temperature of 
~ 1.3 X 10" K seems rather high, however XTE J1118+480 
accretes at rates < 10% I/Edd, while for GRS 1915+105 this 
is ~ I/Edd- Other possibilities for IR contributions include 
optically thick free-free emission (although usually not ob- 
served at such high energies), or optically thin free- free ra- 



diation from an X-ray driven accretion disc wind (Fuchs, 



IMirabel, fc Claret ||2003a| [Castro- Tirado, Geballe, fc Lund 

nggel 



5 DISCUSSION 

Clearly, and perhaps surprisingly, there are large differ- 
ences between datasets 1 and 2, both in the x state of 
GRS 1915+105. While dataset 1 can be flt fairly weU by 
the same model as for BHBs, dataset 2 is marginal, thus 
they do not seem to represent a standard class of spectral be- 
haviour. Interestingly, |Rau fc Greiner] ( |2003l ) also find signif- 
icant variety between the individual x states they observed 
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Figure 6. Multi-wavelength (top) and X-ray band only (bot- 
tom) best fit (Xrod ~1-16) model spectrum (cf. Table Isl SMEDGE 
column), using MNW05-|-Gaussian+smedge. Individualcontribu- 
tions from the MNW05 model are also shown; The light-green 
dashed curve is the pre-shock synchrotron contribution. The dark- 
green dash-dotted curve is the post-shock synchrotron. The pur- 
ple dotted curve represents the thermal contributions from the 
stellar blackbody (below ~ 10~^ keV) and the accretion disc. 



The orange dashed-dotted line above 



10" 



keV represent the 



Compton-upscattered stellar blackbody and accretion disc seed 
photons and the Synchrotron Self-Comptonisation (SSC) of the 
pre-shock synchrotron. The solid grey line is the total MNW05 
model spectrum originating from the jet, the accretion disc and 
the companion, however not forward-folded through the detector 
response matrices and without iron line or reflection contributions 
or absorption due to the interstellar hydrogen column density or 
the smedge model. The red solid lines through the data points 
shows the model flux including all these features. The iron line 
complex, modelled by a Gaussian is shown by the thick green 
curve near 6.4 keV. 



in over 4 years. While the classification done by |Belloni et al.| 
(20001 is based on lightcurves and colour-colour diagrams, 



the sub-indices 1-4 do not refer to a phenomenological clas- 
sification scheme, but only denote their temporal sequence. 
Hence there remains the possibility that the x state still en- 
compasses multiple physical regimes that we may be probing 
here via spectral fitting. 

Alternately, dataset 2 could be in some kind of tran- 
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Figure 7. Multi-wavelength (top) and X-ray band only (bottom) 
best flt (Xred ~l-48) model spectrum (cf. Table Isl REFLECT col- 
umn), using MNW05+Gaussian-(-ref lect. 



sitional state, despite being classified as a typical x state. 
In GRS 1915-1-105 complete state changes can happen on 
very short timescales in comparison to the canonical BHs 
(seconds). Such rapid changes generally occur only between 
states A and B, however. While state C episodes can also last 
for short periods of <100 seconds, the current datasets in- 
volve only the longer (a day or more) intervals. As indicated 
by the radio and NIR lightcurves, it is clearly possible that 
during the dataset 2 observations GRS 1915-H05 was in a 
transitional state and therefore the model, which assumes 
a steady outflow, is no longer applicable. Either way, both 
of the statistically convincing dataset 2 fits (synch and el- 
temp) have other problems in their physical interpretation, 
as discussed above, and it is possible that the model simply 
cannot apply to observations at such high luminosities (see 
Figure [I| . 

For all of the fits presented, two parameters consistently 
stand out by nature of their larger values when compared to 
applications of the MNW05 model to BHBs: the location of 



initial particle acceleration in the jets. 



and the energy 



partition parameter k. As mentioned before, we typically 
find 2acc 1-2 orders of magnitude closer in towards the base 
of the jets in fits of other BHBs in the HS. Similarly, while 
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Figure 8. Multi-wavelength (top) and X-ray band only (bot- 
tom) best fit (x^^^ --^1.44) model spectrum (cf. Tablepl SMEDGE2 
column), also employing MNWOS+Gaussian-fsmedgeTout with a 
reduced magnetic dominance (fc ~ 30 in stead of ~ 550, see table 
i- 



the canonical BHBs often display mild magnetic domination 
over the gas pressure, with k ~ 10, we find at least an or- 
der of magnitude increase is needed to fit GRS 1915-1-105 
plateau states. Finally, the overall powers required (as indi- 
cated by the jet normalisation parameter Nj, see MNW05 
for an explanation) are strikingly larger than the maximum 
observed in other BHBs. Canonical BHBs displaying these 
luminosity levels in I/Edd would long have shut down their 
jets, yet somehow GRS 1915-1-105 seems to be operating on 
a different energy scale. 

These three parameters may well be related. Magnetic 
fields are by now known to play a major role in accelerat- 
ing and coUimating relativistic plasma outflows, and current 
simulations of jet formation favour rather high values of k 
(usually expressed as low values of the plasma /3, denoting 
the ratio to of the magnetic to gas pressure; e.g. |McKin 



ney||2006| [Komissarov et al.|[2007| [McKinney fc Blandford 



2009 1 . It may be that k is positively correlated with the over- 



all power input into the jets. Evidence for such a scaling has 
already been observed in individual fits to broadband data 
of the LLAGN M81* over the course of a yearlong campaign 
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Figure 9. Multi-wavelength (top) and X-ray band only (bottom) 

ll 6kpc column), 
-smedge and em- 



best fit result (Xred ~ 1-01) spectrum (cf. TableKSl 6kpc column). 



for the 1999 data set, using MNW05-|-GaussianH 
ploying a fixed distance of 6 kpc. 



(see, e.g. [Markoff et al.||2008| Figure 22). The overafl mag- 
netic field strength and configuration will certainly infiuence 
the formation of particle acceleration structures, however 
the dependence on total power and other parameters has 
not yet been thoroughly explored. 

A higher level of magnetic domination can also account 
for why the fits favour an electron temperature Te that is a 
factor of two or more lower than the the bottom of the range 
(~ 2 — 7 X 10^" K) found for the canonical microquasars. A 
stronger magnetic field relative to the radiating plasma at 
the base of the jets in our model is synonymous with the 
same being true at the inner edge of the accretion flow, as 
we assume they are directly related. A relatively stronger 
magnetic field would result in relatively more lepton cooling 
in the inner regions, and consequently a lower equilibrium 
temperature. 

In contrast, the geometry of both the base of the jets 
and the cool accretion disc parameters preferred in the fits 
is entirely consistent with the range found in the canonical 
BHBs. The radius of the base of the jets is 3 — 20 rg exactly 
as seen in fits of both BHBs as well as LLAGN. Similarly, 
we obtain best fit values for the inner cool accretion disc 
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Figure 10. Multi-wavelength (top) and X-ray band only (bot- 
tom) best fit result (x^ed "^ ^•^^ spectrum, for the 2005 data set, 



using a reduced low electron temperature Te 
Table [3] eltemp column). 



4 X lO'' K (cf. 



radius rin ~ 2.5 — A.Yrg, which is consistent with such a high 
accretion rate, as well as a spinning black hole as suggested 



by McClintock et al. ( 2006 1. What is not consistent, and dif- 



fers compared to other BHB fits, is that the jet base radius 
in most fits is larger (sometimes only marginally) than rin. 
Since in our model the jets should be launched by the inner, 
hot accretion fiow, this is not internally self-consistent and 
by such criteria, only fits SMEDGe2 and eltemp would sur- 
vive. As the MNW05 model does not yet self-consistently 
solve the jet launching physics from a corona, which is still 
an unsolved problem, we have not made this a hard con- 
straint. 

Although in most BHBs such a small rin would be in- 
dicative of a soft state, this is not true for GRS 1915+105. 
Due to the high accretion rate, a geometry where the ac- 
cretion disc reaches the innermost stable orbit is almost 
expected. The fact that the 6kpc fit has an inner radius 
of less then 1 Vg is most certainly due to the oversimplifi- 
cations mentioned above and similar inner radii were also 
obtained for the canonical BHs (see e.g. MNW05). In fact, 
such small radii appear to be a common occurrence, also 
using other models: analysing 4 years of x state observa- 




Figure 11. Multi- wavelength (top) and X-ray band only (bot- 
tom) best fit result (x^cd ~ 1-15) spectrum, for the 2005 data set, 
using the synchrotron cut-off to fit the X-ray (cf. Table p] SYNCH 
column) . 



tions 
then 



Rau & Greiner ( 2003 1 obtained inner radii of far less 



1 Vg, from the normalisation of the disc contribution 



using the diskbb model. Although this model is known to 
underestimate the inner disc radius by a factor of 1.7-3, due 
to Doppler blurring and gravitational redshift, the obtained 
radii are still unrealistically low. The disc temperatures of 1- 
4 keV they find are higher than ours (0.7-0.9 keV), but again 
the diskbb model is known to produce factor 1.7 overesti- 
mates ( |Lee et aL]|2002| |. Comparing the two plateau states 
to each other, the required disc contribution is much higher 
in dataset 2. For dataset 1, however, disc normalisations 
are high but comparable to those for the canonical BHBs, 
suggesting that, as is true for the HS, the classification of 
individual plateau states should not be based on the lumi- 
nosity. 

In agreement with the findings of many other au thors 



(e.g. Sobolewska & Zycki 



2003 



Rau & Greiner 



20031 who 



have studied the x state, we find that the inverse Compton 
component is always dominant over the direct accretion disc 
contribution in the X-ray regime for our datasets. Although 
in our case the Comptonisation is mainly due to SSC rather 
than purely from thermal accretion disc seed photons. Only 
in one dataset 2 fit (synch) are the X-rays not dominated 
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by the Comptonisation, but the vahdity of this fit may be 
questionable on other grounds as discussed above. 



As brought up as a potential issu e by |Maitra et al.| 
( 2009[) and|Malzac, Belmont, fc Fabian] (|2009[), for high lu- 



minosities conditions at the base of the jet may comprise 
high enough photon densities that pair production can be- 
come important. As these processes are not yet implemented 
in the MNW05 model, we check the pair production and an- 
nihilation rates (hpp and hpa respectively) at the base of the 



jets using the methods described in Maitra et al. (20091, 
and references therein. The results are shown Table |4| If 
the lepton number density n is not larger than iippro/c and 
npaTo/c, where ro/c is roughly the residence time for the 
leptons in the jet base region, we need -hpp < hpa to be able 
to neglect pair production. 

Although our Te ^ 10^"^ is a factor of 5-6 lower than the 
limit deemed problematic for GX339-4, the extreme lumi- 
nosity of GRS 1915-1-105 means that pair production could 
potentially be an issue for some of our fits. The source of 
potential pair production is a high-energy tail above ~ 0.5 
MeV, due to the Comptonisation of thermal accretion disc 
photons. The amount of flux in this tail is directly dependent 
on the normalisation and temperature of the accretion disc 
component, which itself is mainly constrained by the value 
of A^'h and the predicted flux in the tail region. Higher values 
of A^H result in higher accretion disc fluxes to compensate 
for the increased absorption in the soft X-ray band. 

While we fixed the value for Ah in our X-ray fits to 
4.7 X 10^^ cm~'^, an increasing number of works are conclud- 
ing that the column density is actually variable, ranging in 
values from Ah ~ 2 — 16 x lO^'^ cm~^ and potentially linked 



with an intrinsic warm absorber at the highest values (Bel 



loni et al.|2000||Klein-Wolt et al.|2002l|Lee et al.|2002|[Yadw 



2006). Although generally the column density derived from 



X-ray measurements is higher than those obtained using IR 
methods, we find that the highest values of the Ah range are 
not consistent with our model. Values > 7 cause too great a 
decrease below 3 keV, resulting in too much flux in the high- 
energy tail, due to enhanced disc Comptonization. We chose 
the above more "average" value for Ah, in order to be con- 
sistent with p rior X-ray fits of the various x substates (e.g. 
2x10^^ cm~^; Belloni et al. 1997a I and the assumptions of 



[Chaty et al.| ( |1996[ ). For dataset 2 the above value is likely 
too high. Letting Ah free to vary (without correcting the 
NIR dereddening accordingly), while maintaining a higher 
electron temperature of ~ lO'^" K settles on the lower bound 
allowed of 2 x 10^^ cm~^, with a significant improvement in 
Xrcd (from ~30 to ~ 25), mainly due to the resulting re- 
duction in the Comptonised high-energy tail. Therefore our 
results are highly dependent on the value of Ah, and that 
this value is likely different for the two datasets, although for 
the reasons described above we have chosen to use a single 
value. Similarly, the uncertainty in the exact value means 
that too much pair production can be avoided in particu- 
lar by lower values of Ah. Lower values are favoured if we 
allow Ah to vary, and interestingly this is consistent with 



the conclusion by Neilsen & Lee (20091 that the presence 



6 CONCLUSIONS 

Despite the fact that the MNW05 model was originally in- 
tended for application to hard states in canonical BHBs 
only, it appears to well approximate the plateau state in 
GRS 19154-105. However it does not produce convincing 
results in every instance. While some of the parameter val- 
ues obtained are quite extreme, the results for the dataset 
1 appear credible and consistent with what we have found 
in canonical BHBs. Dataset 2 however presents more diffi- 
culties, and a more solid determination of the distance and 
absorption column would go far to help us understand the 
difference between these two plateau states. 

Clearly this work confirms that GRS 1915-1-105 is in a 
much more extreme range of parameter space for an outfiow- 
dominated model, requiring near- or super-Eddington accre- 
tion rates, maximal jet powers and high levels of magnetic 
domination. Although our results confirm previously noted 
plateau state issues, such as the need for a variable Ah, our 
model for the first time incorporates the entire broadband 
and has allowed the comparison between the jet producing 
plateau and hard states. While the baseline geometry seems 
similar, the plateau states of GRS 1915-1-105 are not low- 
luminosity as with HS BHBs, and settle on a range where 
the acceleration of particles occurs much further out in the 
jets, which can be two orders of magnitude further out of 
equipartition in the direction of magnetic domination. While 
these two effects may be linked, the model explored here 
cannot self-consistently address this, but in another work 
we are exploring the links between physical parameters and 
the location of particle acceleration fronts (Polko et al., in 
prep.). A slightly lower electron temperature is also found 
compared to other BHBs, which can be interpreted in the 
context of the higher cooling rates found at GRS 1915-|-105's 
extreme luminosity. 

The main consequence of these differences is that the 
synchrotron component from the outer jets no longer dom- 
inates the soft X-ray band, although a non-negligible X-ray 
synchrotron flux of ~ 10 % the inverse Compton flux below 
~ 50 keV seems required for the statistically favoured flts. In 
comparison, the MNW05 model applied to canonical BHBs 
favours synchrotron emission dominating the flux at least 
up to 10 keV. Interestingly our results are thus qualitatively 
similar to the results found from the blazar sequence (e.g. 
Ghisellini et al. 19981, where higher powers correspond to 



of an absorbing disc wind at high accretion rates should be 
anticorrelated with the presence of jets. 



a decrease in the frequency range where synchrotron power 
peaks. It is clear that time-dependent effects such as cool- 
ing breaks demand further exploration, and we are currently 
working to implement them into more complex models. 

The remarkable differences between two individual 
plateau state observations does raise questions about 
whether the x substates are distinct enough to classify all 
plateau characteristics. While both datasets explored here 
bear all the characteristics of the plateau state, the fits with 
a single model show more variations in free parameters than 
found even between different sources in the HS of canonical 
BHBs. Thus the current classification scheme based solely 
on X-ray colours and timing properties may need to be ex- 
panded based on broadband attributes. 

Our results support a conclusion that, although express- 
ing quite different properties than the HS in canonical BHBs, 
GRS 1915+105 plateau states can still be described by the 
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same broadband model with a steady outflow tied in power 
to the accretion inflow. However the BHB model is clearly 
forced into very extreme ranges, which themselves provide 
some new clues about the relationship between accretion 
rate, jet production and particle acceleration. 

Although GRS 1915+105 is one of the most extensively 
studied BHBs over the last decades since its discovery, we 
are far from understanding this source. At some point in the 
future (assuming the source went into outburst when it was 
first discovered before the end of the century; |Deegan, Com-| 
bet, & Wynn||2009|) GRS 1915+105 wifl invariably retreat 



into quiescence, and should finally yield better insight into 
the connection between its current accretion properties and 
those at sub-Eddington rates. 
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